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ABSTRACT 

It is commonly believed that the magnetic field threading a neutron star provides 
the ultimate mechanism (on top of fluid viscosity) for enforcing long-term corotation 
between the slowly spun down solid crust and the liquid core. We show that this 
argument fails for axisymmetric magnetic fields with closed field lines in the core, the 
commonly used ‘twisted torus’ field being the most prominent example. The failure of 
such magnetic fields to enforce global crust-core corotation leads to the development 
of a persistent spin lag between the core region occupied by the closed field lines and 
the rest of the crust and core. We discuss the repercussions of this spin lag for the 
evolution of the magnetic field, suggesting that, in order for a neutron star to settle to 
a stable state of crust-core corotation, the bulk of the toroidal field component should 
be deposited into the crust soon after the neutron star’s birth. 
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1 INTRODUCTION 


Conventional wisdom states that a strong magnetic field penetrating the core and crust of a neutron star enforces corotation 
between these two components. It is argued that the magnetic coupling timescale is much faster than the viscous Ekman 


coupling timescale and the typical spindown timescale, implying the crust and core spin down in unison. Easson (19791 was 


the first to argue that this issue is far from trivial and attached certain caveats by conjecturing that the degree of crust-core 
coupling is largely determined by the magnetic field topology. As an example, he showed that crust-core corotation on secular 
timescales is not possible along closed lines of a purely poloidal or purely toroidal magnetic field. 


By collating some key previous results (the most relevant being the work of Abney & Epstein (1996l), Melatos (20121 
recently challenged the status quo, arguing that in realistic neutron stars with unmagnetized stratified matter, buoyancy- 
limited Ekman flow can maintain core super-rotation for a timescale of the order of 10^ yr. If true, the consequences for the 
physics of young neutron stars (e.g. Crab-like pulsars and magnetars) are varied and rich. For example, one might expect 
ongoing generation of magnetic field through a persistent dynamo, gravitational wave emission from hydrodynamic instabilities 
and an extra source of pulsar timing noise. 

Spherical, stratified Ekman flow with a magnetic field is currently an unsolved p roblem . [Melatos (20121 argued phe¬ 
nomenologically that stratification lengthens the magnetic coupling timescale (see also [Mendellj l 19981) and concluded that 
for weakly magnetized systems this timescale could exceed the viscous timescale, again allowing for persistent crust-core spin 
lag. Melatos also emphasised the crucial role realistic magnetic field geometry plays in determining the efficiency of magnetic 


coupling (irrespective of field strength), hinging his argument on the aforementioned analysis of Easson (19791. 


It should be emphasized that the two previous arguments, namely, the buoyancy-modified magnetic coupling and Easson’s 
crust-core (no) corotation calculation, are distinct: the former is related to the short-term dynamical response of the core to 
a change in the crust’s rotational frequency, and the latter investigates if an assumed corotating crust-core system is a 
consistent quasi-equilibrium state on secular spindown timescales. This timescale-based distinction is self-consistent, provided 
the short-term magnetic coupling timescale, ts, is much shorter that the electromagnetic spindown timescale, tad = 
and the viscous Ekman coupling timescale, tn. For a neutron star of spin period P and typical internal magnetic field B the 
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magnetic coupling timescale is estimated to be te ~ 5 x 10* G)“®''^(P/10ms) 

which becomes even shorter if one allows for a superconducting core: ts ~ 15 (B/IO^^ G)“®^®(P/10 ms)“^^® s (this is derived 
using the Alfven velocity v% = HcBo/inpp, where He is the lower critical field, and pp is the energy-density of the proton 
fluid). Barring non-superconducting systems with B <C 10^^ G, one always has ts ^ min(tE,tsd), implying the distinction 
between secular and dynamical timescales is valid. 

Investigating the long-term, magnetic field-enforced, crust-core corotation is, therefore, a well defined problem and is the 
subject of this work. Although Easson’s (19791 analytic result is incredibly powerful, it suffers from only being applicable to 


a small subset of possible field configurations, all of which are known to be unstable on dynamical timescales (e.g., Wright 
1973 Braithwaite fc Spruit|2006 Lasky et al. [2011 1. 


In this paper, we fortify and extend the results of Easson (19791 to understand this subtle, yet Important problem of 
magnetic coupling between a neutron star’s crust and core. In doing so, we incorporate linked magnetic fields in which a 
toroidal component weaves through the closed-field line region of the poloidal field. Such ‘twisted-torus’ field configurations 
(e.g., Braithwaite fc Nordlund]2006 Giolfi et al.|2009 l are commonly invoked in the literature as ‘realistic’. 

Our main result Is the following: the existence of any closed poloidal-field-line region in the core of a neutron star eauses 
that region to be magnetically decoupled from the rest of the star. Any region of fluid with open field lines (i.e., those that 
penetrate the crust-core interface) will corotate with the crust, and the two components will therefore spin down in unison. 
A velocity lag will therefore build up between these open-field-line-regions and those with closed poloidal field lines. 

Such a velocity lag between any two fluids is prone to magnetohydrodynamical instabilities that will undoubtably disturb 
the magnetic field equilibrium. We conjecture that such instabilities cause the toroidal component of a twisted-torus-like 
magnetic field to evolve into the crust of the star. The steady-state configuration has the stellar core with only open field 
lines, and hence the entire core corotates with the crust, while the closed poloidal field lines and toroidal field is entirely 
contained within the crust. For newly-born neutron stars with magnetar-like field strengths, evolution to this steady state 
can happen in the first few days to weeks of the star’s life (unless if the birth field is too high, in which case the star spins 
down before the formation of the crust - we quantify and discuss this below). We therefore have a natural mechanism for 
depositing strong toroidal fields in magnetar crusts, thereby supporting theoretical assertions that those strong magnetic field 
components exist (e.g., based on theoretical modelling of the observed temperatures of magnetars ( Kaminker et al.|2006 Ho 
et al.|2012|)). 


The paper is organised as follows. In Section we outline our basic formalism for understanding the quasi-equilibrium 
during spindown and the torque exerted between the crust and the core from the magnetic field. In Section we study 
crust-core corotation for the case of purely poloidal fields. The more realistic case of mixed poloidal-toroidal fields Is discussed 
in Section Section is devoted to the discussion of the astrophysical implications of an incomplete magnetic crust-core 
coupling. In Section we summarize our results and provide pointers for future work. We recommend the ‘fast-track’ reader 
(who is unwilling to navigate the mathematically arduous Sections [^-[^ move directly to Section]^ where they will find the 
astrophysical implications of our results. 


2 CRUST-CORE SPINDOWN 
2.1 Basic formalism 

In this section, we describe the equations that govern the secular dynamics of a neutron star’s fluid core during spindown. 
The solid crust is assumed to rotate rigidly with angular frequency ft{t), which is a slowly decreasing function of time and can 
be identified with the star’s observed spin frequency. Throughout the paper, we work in the instantaneous rest frame of the 
crust, in which the two fluid components of the core, i.e. the neutrons and the proton-electron conglomerate, have velocities 
Vn and Vp respectively. The two velocities are allowed to be distinct in order to account for the likely presence of superfluidity 
in the neutron and proton fluids. However, we shall ignore the likely superconducting nature of the protons and treat the 
magnetic field ‘classically’ in terms of the usual magnetohydrodynamic Lorentz force. 

The long-term dynamics of the charged fluid in the core (protons and electrons) is described by the Euler equatiorQ 

217 X Vp -I- ii X r -I- V^fp = —Fmag - -^Fcpi, (1) 

Pp Pp 

where Fmag is the magnetic force, Fcpi is the coupling force between the neutral and charged fluids, which can include vortex 
mutual friction (see below). Throughout the paper, we use a subscript x = {p, n} to describe proton and neutron fluids 


^ In [Easson] jl97^ , the Euler equation for the charged fluid features an extra force, Fb, attributed to the electromagnetic spindown 
torque ‘ ... exerted directly on the plasma’. We believe this to be erroneous; the electromagnetic spindown torque is exerted on the 
stellar surface (crust) and its effect is mediated to the liquid core through the crust-core boundary and the induced change in B. In other 
words, Fb is already accounted for in F^ag- Strangely, |Easson| | |1979l l drops this force term before concluding his calculation, implying 
our calculations and Easson’s converge. 
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Crust-core spin lag in neutron stars 3 


respectively. In general, the Euler equations feature extra terms related to the entrainment coupling between the neutron and 
proton fluids, however this effect is relatively weak in the core and does not change qualitatively any of the main results of 
this paper. We thus opt for working with a slightly ‘lighter’ formalism with entrainment omitted. It should be noted that the 
omission of entrainment and proton superconductivity takes place only on the macroscopic scale of hydrodynamics while on 
the mesoscopic scale of individual vortices both effects are needed and are implicitly assumed if we are to speak of efficient, 
vortex-mediated, mutual friction coupling between the fluids ( |Alpar, Langer &: Sauls|[l984[ |. 

The left-hand-side of 0 features the inertial Coriolis and Poincare forces and the gradient of 

vl/p = /ip + cj>-l|nxr|", (2) 

which is the combination of the proton-electron chemical potential (per unit mass) and the gravitational-centrifugal potential. 

The neutron superfluid is described by a second Euler equation 

2S7 X Vn -t- X r + VflS'n = — Fcpi, (3) 

Pn 

where Tn is given by an expression similar to 0- 

In the Euler equations we have omitted the inertial acceleration terms, which are negligibly small when Vp and Vn are 
small. This approximation implies a nearly corotating crust-charged fluid system. The state of exact corotation is described 
by Vp = 0, in which case the proton and neutron fluids are rigidly rotating about the same axis as that of the crust (which we 
identify as the 2 :-axis) with angular frequencies as measured in the inertial frame fin = flp = fl- Then, for the nearly-corotating 
crust-core system we can assume flp, fin ~ fl without fin and flp necessarily being equal. 

Equation 0 describes the departure of the system from a ‘background’, time-independent state of exact corotation, 
where this departure is driven by ji. As a consequence, all quantities appearing in 0 are perturbations with respect to 
the corotating state. For instance, the magnetic field can be decomposed into a fixed background part, Bq, plus a small, 
spindown-induced perturbation, b = dB, i.e. 

B « Bo + b. (4) 


The magnetic force can therefore be expressed to leading order 


-V 


Bob 

Att 


-b — [(b-V)Bo-b(Bo- 

47r 


V)b]. 


( 5 ) 


The nature of the coupling force, Fcpi, is determined by the interaction between the superfluid’s quantized vortices and 
the other components of the star (e.g. protons, electrons, the magnetic field). The simplest form for this force is provided by 
the Hall & Vinen mutual friction force ( Hall fc Vinen|1956a|b| | 

Fcpl = 2flnPn { S [z X (z X w)] -b S'(z X w) } , (6) 

where w = Vn — Vp is the velocity lag between the neutrons and the charged particles, B and B' are mutual friction coefficients 
(which we take to be uniform throughout the star) and a ‘hat’ denotes a unit vector. This form of the force assumes a vortex 
array aligned with the common spin axis. According to our setup of the problem we have, 

w « w(fln - flp)(p = urQ.npip, (7) 

where zu = r sin 0 is the usual cylindrical radius. The mutual friction coupling force becomes 

Fcpl = —2ti7pnflnflnp (Hp + H'to ) . (8) 

The most commonly considered mutual friction mechanism in neutron star cores is the scattering of electrons by the magnetic 

Andersson et al.[2006(. 


field of individual vortices, in which case B' B ^ 10 ^ (Alpar, Langer & Sauls 


1984 


It should be noted that if (for whatever reason) superfluidity is absent, Fcpi would represent the force due to inter-particle 
collisions and the core would effectively behave as a single-fluid system with one common velocity, described by eqn. (Q with 
the appropriate readjustment of the density and pressure. The entire analysis of Sections [2||^ still hold in the absence of 
superfluidity. 

Following Easson (19791, it is easy to show that the Coriolis force can also be omitted in eqns. 0 and 0 - The argument 
goes as follows: the Poincare force in 0 is the term driving the magnetic field evolution. We estimate the strength of the 
magnetic field perturbation, 6, by balancing the magnetic and Poincare forces: 

" P 


PpHP 


Bo 271 


VA ) 


tsd 


( 9 ) 


where Xp = Pp/p is the proton fraction (p = p„ -b Pp is the total density), t^o = H/|H| is the characteristic spindown timescale, 
v\ = BoCtzp is the Alfven speed, P is the spin period and R is the stellar radius. 

Our perturbative approach is justified as long as b Bq or, equivalently, when 

1/2 /I 1/2 


Bo » Pc: 


10 " 


10 ms 


lyr 

tsd 


( 10 ) 
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4 Glampedakis & Lasky 


where we have assumed canonical values for the stellar parameters, R — 10® cm, p = gcm“®, a;p = 0.05, and we have 
normalised P and tad with new-born, rapidly spinning down systems in mind. This condition is indeed expected to be satisfied 
by the vast majority of neutron stars. Had we used the superconducting expression for the Alfven speed, V\ = HcBo/inpp 
(where He is the so-called lower critical field), we would have obtained a negligibly small threshold for Bcrit ( |Easson|1979[ )). 

We can now show that 6 <C Bo implies a negligibly small Coriolis force as compared to the Poincare force. From the 
perturbed induction equation for the magnetic field we have. 


The ratio of the two forces Is then 


atb = V X (vp 

xBo) 

-DO Csd 

(11) 

2Sli!p 

QR 

Vptsd b 

R ^ Bo ’ 

(12) 


Implying we can safely omit the Coriolis force in the Enler equations 0 and 0. 


2.2 Easson’s argument against corotation 

If we assume the magnetic field is symmetric with respect to the spin axis, both the background and perturbation, Bq and b, 
can be decomposed into poloidal and toroidal components 

Bo = Bp + h = hv + b‘^<f>. (13) 

The azimuthal component of 0 can be combined with equations 0 and 0 to obtain 

B-p • V {vjb‘^) + bp • V (zuBq) = 47rW^ ^ —pptl + 2Bpn^n^np ^ . (14) 

which, for the case of a purely poloidal background field, Bq = Bp, becomes 

Bp • V^ —Pptl -I- 2Spnflnllnp ^ . (15) 


Easson (19791 argued against crnst-core corotation based on equation (15l. Snppose the field lines of Bp form a closed loop 


somewhere in the core. The integral of the left-hand-side of (151 along that loop vanishes (assuming b^ is not multivalued) 


while the integral of the right-hand-side is finite. Therefore, equation ( |15| does not lead to a physically acceptable solution 
for b, signalling the failure of the corotation approximation along closed poloidal lines. 

The assumption of a purely poloidal background field is overly restrictive, especially given that such fields are unstable 
on Alfven timescales (e.g., Wright 1973 Braithwaite & Spruit 2006| Lasky et al. 20111, which is much shorter than the 
secnlar timescales we have considered heretofore. Mixed fields that include poloidal and toroidal components, where the 
toroidal component threads the closed-field-line region of the poloidal field in the stellar core, provide dynamically stable 
configurations (e.g., Braithwaite & Nordlund 2006 Akgiin et al. 20130 Although such ‘twisted-torus’ configurations are 


much studied, their effect on the rotational dynamics of a neutron star core-crust system are not known. In this paper, we 
calculate the effect various magnetic field confignrations have on neutron star crnst-core coupling, which relies on calculating 
the torque between these two components — a task to which we now turn. 


2.3 The crust-core torque 

The magnetic coupling between the crust and core implies a non-zero traction (force per unit area), t, at the base of the crust. 
The differential torque per unit area exerted on the crust by the core is 

M — r X t. (16) 

The azimuthal component of the traction is responsible for a non-vanishing torque along the spin axis. This is built from the 
azimuthal and radial magnetic field components. 


The corresponding torque per unit area is 


tv = —B‘^BP 
47r 


A7" = —B'^BP 

47r 


(17) 


(18) 


It shonld be noted that the torque is calculated at the base of the crust, at radius r — Rc, in which case zj = zjc = Rc sin0. 


In the problem at hand, the magnetic field is perturbed as a response to the spindown. Equation (181 can be expressed 

J^^ = ^{b^B^p+PB^), (19) 

47r 


^ [Lander &: Jone^ | |2012[ | and [Mitchell et al.| | |2015| | have shown that such fields are unstable in non-superfluid barotropic stars, a result 
that is not applicable here. 
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Crust-core spin lag in neutron stars 5 


where we assume that the background magnetic held can have both poloidal and toroidal components at the crust-core 
boundary. 

The total torque exerted on the crust, A^cc, is given by the surface integral of A/”^ over the crust-core boundary (which, 
to a good approximation, can be assumed to be spherical): 


iVcc = / dSAf^ = ^ 
47r 


L 


[ dSwcib‘^Btt+b^B^). 

J Rc. 


( 20 ) 


The solution of the Euler equation provides b'^ and 6’’ in terms of Q. and fix which, upon substitution in (201, eventually leads 
to an equation of the form 


Ncc — —dCl A^cpi- (21) 

where Ncpi denotes the mutual friction coupling term and the coefficient / plays the role of a moment of inertia (see below). 
The torque, A^cc, enters the equation of motion for the entire solid crust: 


/cr^i — A^cc + A^em, 


( 22 ) 


where is the crustal moment of inertia and Nem 
and rearranging gives 


is the spindown torque due to the exterior magnetic field. Inserting (211 


(Icr + 7)12 — A^em + A^cpl. (23) 

Equation ( |23[ ) can be viewed as the equation of motion of the combined system comprising the crust and the fraction of the 
core with moment of inertia I. In other words, I represents the fraction of the charged fluid in the core that is efficiently 
coupled to the crust and therefore spins down in unison with it on a secular timescale. 

In general, I may not coincide with the total moment of inertia, 7p, of the proton-electron core fluid. When I ^ Ip, the 
magnetic field cannot enforce corotation between the crust and the entire charged core. In that case, a portion of the core 
fluid will fail to follow the spin evolution of the crust and, barring the action of some other non-magnetic crust-core coupling 
mechanism, a spin lag will develop between the uncoupled fluid and the rest of star. The rest of this paper is concerned with 
calculating I for different magnetic field conhgurations, and hence determining which portions of the stellar core corotate with 
the crust. 


2.4 A note on symmetries 


Before moving on to the main part of the paper, i.e. the detailed calculation of the crust-core torque, Afcc, and the moment 
of inertia, 7, we need to discuss the symmetry of the magnetic perturbation, b, with respect to the equatorial plane. Given 
the overall axisymmetry of the system, the components of b can be either symmetric or anti-symmetric under reflection with 
respect to the equatorial plane, i.e., 2 —>■ — 2 . An inspection of equation ( |14| for the core spindown quasi-equilibrium reveals 
that the right-hand-side term is reflection-symmetric. Consequently, it only couples to those components of b that would 


make the left-hand-side terms of (14l equatorially symmetric when combined with a given background held Bq. 


The magnetic helds considered in this work are characterized by an anti-symmetric radial component, Bq{z) = —Bq(—z), 
and symmetric Bq(z) = Bq(—z) and Bq{z) = Bq{—z) components. This symmetry property of Bo implies that only an 
antisymmetric b‘^ and a symmetric 6’’ can make a contribution to the integrated torque Ncc (see equation (|19[)). It also means 
that, in order to calculate the torque, it suffices to integrate over just one hemisphere (e.g., the upper one). 


3 CRUST-CORE COROTATION: POLOIDAL FIELDS 

We shall hrst study the problem of crust-core corotation for the case of a purely poloidal background Held, i.e. Bo = B-p, 


without closed lines anywhere in the core. The relevant equation for the spindown of the core is eqn. (151: 

Bp • V(t*76'^) = ^ —Pptl -I- 2SpnS2nI2np j . 

Quite generally, the background Held can be written in terms of an axisymmetric stream function, i/)(r, 6), as 

Bp = Ntj) X V<p = —Nip X (p. 


(24) 


(25) 


By construction, Bp is tangent to the ip = constant surfaces (the introduction of the stream function also makes V • Bo = 0 a 
trivial identity). Hereafter we shall only consider the special case of a general dipole Held. The corresponding stream function 


Ip = Af{r) sin^ 0, 


(26) 
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6 Glampedakis & Lasky 


where ^ is a constant amplitude and /(r) is an arbitrary fnnction (modulo the requirement of regularity at the origin). The 
held components can then be expressed in spherical coordinates, 


Bq = cos 9, 


Bq = -A^^sine, 


where a prime denotes differentiation with respect to the argument. 

/' 


Inserting these in eqn. (241, 


2^ cosOdrT — ^ sin^^eT = ^ ( ~Pp^ + 2.BpnIInIInp j 


(27) 


(28) 


where we have introduced the toroidal function 


T{r,e) = vjy<’. (29) 

The resulting partial differential equation is not, in general, amenable to a solution through separation of variables. In |Easson| 
( |1979[ ) this equation is solved for /(r) = 1/r combined with the extra simplihcation of a uniform density star. This choice 
of magnetic field would be a suitable one for the exterior space but is clearly an unphysical model for the field in the stellar 
interior. An example of a separable field that is also regular at r = 0 is /(r) oc however this also requires a uniform density 
profile and therefore is of limited physical interesl[^ 


3.1 The simple case of a uniform field 

The only separable solution admitted by realistic, non-uniform density neutron star models is /(r) oc r^, corresponding to 
the uniform field Bq = Bqz with Bq = 2A and BJ = 2Acos6', Bg = —2Acos^. However, the cylindrical symmetry of the 
field suggests that it would be easier to study the problem in cylindrical rather than spherical coordinates. Indeed, the Euler 


equation (241 in cylindrical coordinates becomes 

Bodz(pjb'^) = 47rW^ ^ —pptl + 2Bpnllnflnp j . 


(30) 


This can be integrated without any difficulty (the integration constant is fixed by the requirement of an equatorially anti¬ 
symmetric b'^), 


b‘^{zu,z) = fi J dz Pp + 2 Bi}n^np J dz' Pn'^ ■ 


The total crust-core torque from eqn. (201 is given by 


Ncc = 2nR, 


f 

Jo 


dvj- 


Blb^ 


(B2-tu2)l/2 4 tp 

and the radial component of the background field is (in the upper hemisphere) 


Bq = Bq cose = Bo^ = B( 




p p ’ 

Jtc -ttc 

where Zz{vj) = y/R^ — is the «-coordinate of a point on the spherical crust-core boundary. 
Given these expressions, the torque becomes 


A^cc = 47r 



" . c 

/ dzuuj'^ ( 

-n / 

Jo ' 

V Jo 


dzpp + 2Bflnfln 


, J dzpn ^ 


= -IpQ + 2BRQ.A 


where 


7x = 47r 


dzuzu [ 

Jo Jo 


Zc{-^) 


dz 


(31) 

(32) 

(33) 

(34) 

(35) 


is the total moment of inertia of the x = {n, p} fluid in the core. The equation of motion for the crust, eqn. (231, takes the 
following rigid-body form: 


( 7 cr + 7 p)f 2 — A^em + 2 BIn^n^n 


(36) 


The comparison of eqns. (231 and (361 allows us to conclude that the uniform dipole magnetic field, Bq = Bqz, leads to 
I — Ip, which means that the entire proton-electron fluid of the core is coupled and corotates with the crust. This result 
was first obtained by Easson (19791. We can show that the same property is true for the earlier case of the non-uniform field 
f{r) = in uniform density stars — see Appendix [a| for a detailed discussion. 


® It should be pointed out that the radial profiles / oc (r^,r^) represent the solutions of the dipolar Grad-Shafranov equation that 
describe the hydromagnetic equilibrium in uniform density neutron stars ( Ferraro|1954^. 
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Crust-core spin lag in neutron stars 7 


3.2 Solving the corotation problem in magnetic coordinates 


The ease with which the uniform field case was solved was a consequence of using a coordinate system adapted to the geometry 
of the field lines. This motivates us to use magnetic coordinates for a general, axisymmetric poloidal held without closed lines, 
with the expectation of a troubleless solution for the perturbed toroidal magnetic held Moreover, the assumption of 
uniform matter is hereafter permanently abandoned, leaving the stellar equation of state completely general. 

The magnetic coordinates comprise the stream function, ip (any other function T’(^) could be equally well used), the 
length along a held line, y, and the previously used azimuthal angle, ip. The set {'ip, x, p} dehne an orthogonal curvilinear 
coordinate system with a line element 

dl^ — h^dip'^ + C^dx + hTpdp^, (37) 

where the metric functions ,/i,,, are all (^-independent. Appendix [B| provides the necessary mathematical supplement 

on magnetic coordinates. 

Expressed in the new coordinates, eqn. ( |15| l for the spindown of the core charged huids becomes, 

So, 


-d^T = ^ —Ppfi -f 2f3pnflnflnp ^ ■ 

As anticipated, this can be trivially integrated; 

Ti'ipiX) = / dx -\-S7rBQMnp / 

Jo -DO Jo DO 

where we have taken x = 0 at the equatorial plane and enforced an anti-symmetric (see Section 2.41. 


(38) 


(39) 


For the remainder of this calculation, we need two geometric identities concerning our magnetic coordinates, the calcu¬ 
lations for these are presented in Appendix[^ Firstly, the surface element, dS, of a spherical surface in magnetic coordinates 
is 

1/2 


dS — 


1-b 


ttp Sc 

hip Re 


dpdpi, 


(40) 


where Xc = Xc{Rc,'ip) is the x-coordinate of a point on the sphere with radius r = R^. Secondly, for a general dipole held, 
'ip = Af{r) sin^ 9, we have the exact identity 


C-x) 


^^ f hxdpr 

yhipd^r 


1 1/2 


= 1 . 


The crust-core torque from eqn. (|20|) can now be expressed as 

dpdip 


iVee = -r'- f dSvoeb'^BE = ^ f dpdiP 1 -b ( 

471 47ry^^ \hipdxRcJ 


1 1/2 


{hiphipB^T)^ 


Using Bq = Bc(f ■ x)xc> together with our earlier solution for T in eqn. (391, Nee becomes 

1/2 




[h.f,hipBo{f ■ x)]xo i - 




K 


h It 

dx-^PpW^ + 2Z3nnf2np / dx'-^pn-!^^ 

t >0 Jo 


(41) 


(42) 


(43) 


The integral can be simplihed using Bq = l/{hiphip) (see eqn. (B4l), and the torque formula can be expressed in the desired 
compact form 

Nee = —fpil -\- 2Zlflnf2np7n, (44) 

where the effective proton and neutron moments of inertia are dehned as 




1/2 


C ■ x)k 


L 


XcW 


dxhxhy^hppx'oo . 


(45) 


where x = {p, n} for protons and neutrons respectively. On the other hand, the x-huid’s total moment of inertia in magnetic 
coordinates is given by 


/x = 


/ dxdShx{r-x)p^m'^ 

Jr<Rr~ 


which, after carrying out the (/)-integration and using the upper-lower hemisphere symmetry, become^ 


7x = 471 




XciP) 


dxh,phiphx(r ■ x) 


1-b 


/ hxdp: 
V hpdx 


1/2 


2 


(46) 


(47) 


It is worth noting that, in the limiting case of cylindrical coordinates such that we have = 1 and 

r ■ x = cosO = zjr., and |[47| reduces to (|35|l. 
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Figure 1. Schematic configurations of mixed poloidal-toroidal magnetic fields. The outer crust and crust-core interface are the red and 
blue circles labelled R and Rc respectively. Left panel: closed field lines in the core — the thick blue field line, labelled A, is the outer-most 
poloidal field line that is closed in the core, and the red magnetic field line, labelled B, is the outer-most closed field line in the entire star. 
For the twisted-torus model, toroidal field is confined to the region enclosed by B. The black dots represent the neutral line (the point 
where the poloidal magnetic field component is zero). Right panel: no closed field lines in the core - the only closed field lines penetrate 
the crust-core interface, with B being the outer-most closed poloidal field line. Again, the twisted-torus model has toroidal field confined 
to the region enclosed by B. 


It follows from (451 and (471, together with the geometric identity (411, that 

/x - /x. 


(48) 


We therefore conclude that a general axisymmetric dipole poloidal field (with no closed lines) couples the entire charged fluid 
core to the crust, thus enforcing corotation between these two components. 

This result does not support Easson’s claim that poloidal fields (with open lines) do not always enforce crust-core 
corotation (Easson 19791. It is not difficult to find the reason behind this disagreement: as we have already pointed out, the 
counter-example calculation provided by Easson makes use of the unphysical interior dipole field / = 1/r. 


4 CRUST-CORE COROTATION: MIXED POLOIDAL-TOROIDAL FIELDS 


The preceding purely poloidal background field model is an ideal testbed for understanding and solving the magnetic crust-core 
corotation problem. However, based on arguments about dynamical stability, it is commonly believed that realistic neutron 
stars are likely to harbour a field of mixed poloidal-toroidal geometry. Therefore, it is highly desirable to extend the previous 
analysis to encompass the case of a mixed magnetic field. 

The assumed geometry of the field is displayed in Fig. [^and makes contact with the twisted torus configuration discussed 
earlier. More specifically, the poloidal field threads the entire star and has a neutral point (fine) somewhere in the core. The 
toroidal field is confined inside the star, occupying the region around the neutral point, and is bounded by the last closed field 
fine inside the star, i.e., fine B in the figure. On the other hand, the fine A marks the last closed field line in the core. As a 
result, the mixed poloidal-toroidal field fines that thread both the crust and core are the ones between the lines A and B. 


As we have already seen, the spindown quasi-equilibrium for a mixed background field, Bq, is described by eqn. (141, 

(49) 


Bj 


V{voh‘^) + b-p • V{vjBq) = 47rw^ ( —Pptl + 2I?PnHnfInp ) 


Since bp is axisymmetric, we can introduce a new stream function l3{r,9): 

bp = V/3 X V<p. 

Expressed in magnetic coordinates, where /3 = P{ip,x), 


bp = 


h{p hip 


hip h^ 




We also introduce a stream function X for the background toroidal field, 

= TV VP. 


(50) 


(51) 


(52) 


© 0000 RAS, MNRAS 000, 000-000 










































Crust-core spin lag in neutron stars 9 

It is well known that for an axisymmetric magnetic equilibrium the two stream functions are related as X — XCi) (^-g- 


Glampedakis et al. (2012|). 


In terms of these new functions, we can write (491 as (recall that B-p = 1/h^hip, see Appendix E 

d^T — X'd^P = ^ —Ppfi + ^ • 


(53) 


It is straightforward to see that (531 breaks down when integrated along a closed field line in the core. Along such a line X' 


is constant and therefore the integration leads to 

0 ^ ^ d^zu ^ ppH A 2ZIpnHnHnp ^ , (34) 

which is a contradiction. This verifies the assertion discussed earlier in Section [2^ that closed field lines are incompatible with 
the assumption of crust-core corotation. 

Next, we consider ( |53[ ) along an open field line. Its integration yields 

rx rx 

M(x, V') = r(x, V’) — 2i'/3(x, V*) = —47rn / dx'h^ih^h^ppW^-\-8nBtlnptdn dx'h^/hph^ppU^'^, (55) 

Jo Jo 

where we have used the assumed equatorial symmetry and anti-symmetry for b'" and 6'^, respective!}]^ 

For the crust-core torque we have (remember that this is the contribution of the upper hemisphere multiplied by two): 


1 


dS 


= - 


{M + dCXP)}^^^-X^{f-^) 


(56) 


In taking the surface integral, we need to make sure that the range of ip is limited to the open lines across the crust-core 
interface. That is, ip should ‘run’ between the value ipz of the symmetry axis and ip a of the line A, see Fig. At the same 
time, the range for the ^-integration of the terms that depend on the toroidal field X is the subset ips < ip < ip a of the 
previous range, with ipB corresponding to the boundary line B of the toroidal field region. 

It is straightforward to obtain, 


Ncc — — Ip^ A 2ZIIInIInp/n A 

where we have defined 




1 1/2 


d^{XP){f ■ x) -XdxP^if ■ Ip) 


/x = 47r 


j-'^A 

/ ' 

J 


dip 


1 A 


f t^x tJ'P 

y /It/ 


1 1/2 


[r-Xix. 

Jo 


Xc(V>) 


dxhxhph^p^vj . 


(57) 


(58) 


This can be compared with the moment of inertia appearing in the earlier calculation of the poloidal field with open lines. 
The difference in the integration limits implies that /x < 7x, with the equality holding when there are no closed lines in the 
core or, equivalently, when the poloidal’s field neutral point is pushed outside into the crust. 

So far the discussion in this section applies to the case of an arbitrary twisted torus-like held, but in order to proceed 
further we make the assumption of a dipole poloidal held Bv, i.e. 


Ip = Af(f) sin 9, 


i?’ 


(59) 


where for convenience we have taken / to be a dimensionless function. In Appendix]^ we obtain the following results for the 
geometric factors appearing in (571 and (581: 

-li/2 


1 A 


Be 

h-ipdxBc 


1 1/2 


(7-x) = l, ^ 


1 + 


Rc 

hijjO^Rc 


(f • -0) = ± 


2Af cos 9 ’ 


(60) 


where the explicit form of the function G{r) is given in Appendix]^ The sign arbitrariness of the second equation comes from 
the fact that the inner product r ■ ip can switch sign along a given held line (open or closed). 


We then have from (571 and (581, 


L = 471 


/ dip / 

Jijj^ Jo 


XaW 


dxhxhptlip p^taj a 


and 


Nec — —7pf2 A 2B^nt^nptn A X{ipA)pJlXc^ V’a) a 


Gc 

2Afe 


rpA 

/ ' 


dipx 


idxlJ)x. 
cos 6 


( 61 ) 


(62) 


where Gc,/c are the values of G, / at r = Re- We have also used the fact that X{ipB) ~ 0 for the assumed magnetic held 
geometry. With the help of (see Appendix [b|) 


X = G cos 9 


XeilP) = Gc ( 1 - ^ 


1/2 


(63) 


® A symmetric b’’ entails an antisymmetric function this becomes evident since in the limit 2—^0 we have b’f —> fe’’ and dx —t 5z. 
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10 Glampedakis & Lasky 

we can obtain the equivalent expression, 


N^c = -Iptl + 2B^ln^npIn + I{'tpA)P{Xc, '^/’a) T TVjV / 

2Ajc 


1 - 


jtp_ 

AU 


- 1/2 

AdxP)x.- 


(64) 


In order to proceed further, an independent solution of the form /3 = aifl + a 2 flnpfln must be provided; such a result can 
only be obtained by solving the remaining poloidal components of the Euler equation 0- These equations can be found in 
Appendixhowever, solving them analytically does not appear to be possible. 

This difficulty did not crop up in our earlier poloidal field analysis because in that case the (p-Euler equation, the torque 
Ncc and 6“^ itself decouple from bp and the poloidal Euler equations. In practise, this means that a purely poloidal Bq field 
Is compatible with the toroidal spindown ‘mode’ (6‘^, bp, dTp) = (6'^, 0,0). 

A purely toroidal spindown mode is no longer admissible when Bp yf 0 but, nevertheless, it is still possible to find a 
mixed toroidal-poloidal mode that would allow a complete calculation of the crust-core torque The desired b is simply 
the superpositior0 

^ = Kr + ble + bl(p, (65) 


where ‘a’ and ‘s’ stand, respectively, for antisymmetric and symmetric. 

As a result of the imposed symmetries we have that (i) the ~ Bq term in the crust-core torque (201 integrates to zero 
and (ii) the second left-hand-side term in (141 is anti-symmetric (as opposed to the other symmetric terms) and therefore 
should be zero, l.e. bp • \/{wBq) = 0 —>■ = 0. In other words, this spindown mode is characterized by a poloidal 

perturbation that points along the background field lines, bp = b^X- At the same time the perturbation (651 also leads to a 
pair of non-trivial poloidal Euler equations for /? and ^'p. 

Adopting (651 as a spindown ‘mode’ implies that we can remove all /?-dependence from the general result (621 (or (64l), 
leaving 

iVcc = -Iptl + 2Z3BnfInp/n, (66) 

as the main result for the total crust-core torque for the case of a mixed magnetic field. According to this, the magnetic 
field couples the crust with exactly that portion of the core that is threaded by the open magnetic field lines. Assuming that 
poloidal field has a neutral point somewhere in the core (as in the left panel of Fig.[^, the toroidal-shaped region between the 
field lines A and B (i.e. the shaded region in Fig. remains magnetically decoupled and of course makes no contribution to 
the torque. The same property can be seen to hold true for the more general case represented by eqn. (64l. The astrophysical 
Implications of this remarkable result are discussed in the following section. 

Finally, it is worth remarking on a fundamental difference between the above argument regarding the crust-core torque 
and Easson’s no-crust-core corotation argument. We remind the reader that Easson (19791 showed that a crust-core coro¬ 
tation equilibrium cannot exist along closed poloidal field lines in the core, as the governing equations lead necessarily to a 
mathematical Inconsistency when corotation is assumed (see Section 2.21. Our analysis, facilitated by the use of magnetic 
coordinates, has shown that the closed field lines of poloidal/mixed magnetic fields do not couple at all with the crust and as 
a result the crust cannot exert any torque on the region occupied by the closed lines. These two anti-corotation arguments are 
clearly complementary; our analysis explains that closed field lines can violate the corotation assumption as they are decoupled 
from the crust. This subtle issue may be purely academic; as we discuss below, we expect the spin-lag generated between the 
open- and closed-field-line regions to render the latter region dynamically unstable, implying crust-core corotation can only 
be established in the absence of closed field lines in the core. 


5 FAILURE OF LONG-TERM CRUST-CORE COROTATION: IMPLICATIONS 
5.1 Incomplete crust-core coupling: how likely is it? 


The results obtained in the preceding sections strongly suggest that the extent to which the liquid core can follow the secular 
spindown of the solid crust is largely decided by the geometry of the magnetic field in the core, in agreement with the original 
conjecture put forward by Easson (19791. A key question raised in that early work concerned the ‘size’ of the class of magnetic 
fields that are unable to enforce a full crust-core coupling, and hence allow for the build up of a persistent spin-lag between 
part of the core and the rest of the star. 

This paper provides a definitive answer to that crucial question. As we have seen, an arbitrary dipole poloidal field with no 
closed field lines in the core does enforce complete crust-core corotation, in agreement with ‘conventional wisdom’. However, 
the presence of purely poloidal magnetic fields in neutron stars enjoys very little theoretical support since these fields are 


® In terms of a multipole expansion in Legendre polynomials Pi{cos 9) this superposition would correspond to an £ = 2, 4, 6,... sum for 
b'^ and an £ = 1, 3, 5,... sum for b-p. 
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inherently unstable (e.g. Wright (19731; Flowers &i Ruderman (19771; Braithwaite (20071; Lasky et al. ( 2011[ )) (the same is 
true for purely toroidal fields ( Thyler|l973 Braithwaite|2006 1). Instead, hybrid poloidal-toroidal ‘twisted torus’ fields are the 
ones commonly attributed to realistic neutron stars and the ones that have attracted most theoretical work in the recent 


years (e.g. 

Braithwaite & Spruit ( 

20061; Colaiuda et al. (20081; Ciolfi et al. (2009|); Mastrano et al. (20111; Glampedakis et 

al. (20121; 

Lander & Jones (|2012 

1; Ciolfi & Rezzolla (20131; Fujisawa & Kisaka (|2014|)). Such fields are exemplified by the 


crust for sufficiently strong toroidal components ( Ciolfi et al.||2009 l) and, as we have shown in Section]^ the part of the core 
occupied by these lines remains magnetically decoupled from the crust (and the rest of the core). 

This result implies that the aforementioned class of magnetic fields is ‘large’, in the sense that such fields should be 
considered likely to occur in real systems. These include the most commonly used twisted torus models, but not necessarily 
fields for which the magnetic field is misaligned with the rotation axis (see section]^. In what follows we discuss the spindown 
evolution of a neutron star endowed with a twisted torus-like magnetic field and speculate on the implications of our results. 


5.2 Implications for the magnetic field evolution 


To what extent the crust and the core can be interlocked in long-term corotation is most relevant soon after a neutron star’s 
birth. Assuming a standard cooling scenario, the crust is expected to form approximately a day after the neutron star is 
formed, when the temperature has dropped to T ~ (few) x 10® K ( Shapiro fc Teukolsky|[l983 1. During the same period, the 
star retains most of its initial kinetic energy (this is true even for newly-born magnetars unless their surface dipole field 
exceeds approximately G) and therefore most of the subsequent electromagnetic spindown will take place when 

the star has a well defined solid crust and a liquid core. The onset of neutron superfluidity in the core is likely to ensue much 
later (~ lyr— 100 yr), when the star has cooled below T ~ 10® K— (few) x 10® K, respectively (Shternin et al. 


2011 


Page et al. 


20111. On the other hand, the formation of the superconducting proton condensate is expected to take place much earlier, at 
T > 10® K ( Shapiro fc Teukolsky|[l9^ K 

Once the crust is well formed and has begun to be torqued by the anchored magnetic field, two mechanisms come into 
play in an attempt to enforce crust-core corotation: the magnetic field and hydrodynamical viscous Ekman pumping. The 


latter mechanism has been extensively discussed in the neutron star literature, most notably in Easson (19791; Abney & 


Epstein (19961; Melatos (20121, see also Greenspan (19681. As recently emphasized in Melatos (20121, the combined effect of 


stratification and compressibility in neutron star matter severely inhibits the efficiency of the Ekman mechanism, resulting in 
corotation timescales much longer than initially anticipated ( Abney fc Epstein|1996[ ). For typical neutron stars parameters, a 
viscous coupling timescale as long as ~ 10® yr is expected ( Melatos||2012 K The magnetic coupling mechanism is much faster 
than this (see also our discussion in Section [^; we can therefore frame our discussion without having to consider viscosity, 
consistent with the model presented in this paper. 

A day-old neutron star is likely to start its secular spindown evolution being in stable hydromagnetic equilibrium with a 
magnetic field geometry that resembles the ones shown in Fig. Gonsider first the case of the right panel of Fig. In that 
scenario, the poloidal neutral point is located in the crust, as is the bulk of the toroidal field. The entire core is threaded by 
open (poloidal and toroidal) field lines and, as we have shown in this paper, the magnetic coupling can enforce corotation 
between the crust and the total amount of charged fluids in the core. Our analysis and results apply equally well to a system 
with normal or superfluid neutrons, in the former case the implication being that the magnetic field directly couples all the 
fluids to the crust, while in the latter case vortex mutual friction is very efficient in enforcing corotation between the superfluid 
and the protons. 

More interesting (and perhaps more probable) is the case of neutron stars ‘starting’ with the magnetic field equilibrium of 
the left panel of Fig.[^ where the poloidal neutral point and the bulk of the toroidal field are located in the outer core. In such 
a case, the toroidal-shaped region bounded by the closed line A (depicted as a shaded region in Fig. remains magnetically 
decoupled to the crust and, as the rest of the star spins down, a spin lag will begin to build up, with the toroidal region 
retaining its initial ‘fossil’ spin frequency. 

It is unlikely that this state of persistent crust-core spin lag can remain dynamically stable for long. The velocity jump 
in the proton fluid across boundary A will give rise to an increasingly strong azimuthal current sheet A Jip and an associated 
poloidal Lorentz force AJ(f> x 'Q-pjc. We would expect |AJ/J| ~ Afl/fl where AQ{t) ~ Qt is the spin lag. This means that 
the induced Lorentz force can seriously disturb the background equilibrium after a time t ~ tad = D/|f2|. 

The physical conditions along the same boundary field line A are also suitable for the onset of the classical Kelvin- 
Helmholtz instability. This should take place once the spin lag exceeds the threshold ADkh ~ va/R, where va is the average 
Alfven speed of the local toroidal field (see Ghandrasekhar (19811, Chapter XI, p. 507). The corresponding timescale for 


the onset of the instability, Ikh, can be estimated to be (here we use canonical values M = 1.4Mq and R = 10km for the 
stellar mass and radius): 


tKH 

fsd 


VA 

nR 


5 X 10 


-2 j 




\ 10^® G/ \ 10ms 


(67) 
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12 Glampedakis & Lasky 


suggesting that it can set in much sooner than the lapse of one spindown timescale. The inclusion of proton superconductivity 
can significantly raise the above ratio, va/Q.R ~ 0.2 G)^^^(P/10ms), but this is still well below unity provided 

P < 10 ms, assuming a magnetar-strength field Bq . 

We thus have good reason to expect that the initial persistent crust-core spin lag will lead to some kind of magnetohy¬ 


drodynamic instability after a time (assuming a standard spindown torque, see e.g. Shapiro & Teukolsky (19831) 

-2 


t<tsd~ 4.7 


Bd 


10 ms 


( 68 ) 


^lOiSG, 

where Bd is the dipole surface field. Among all systems, newly born magnetars appear to be the most favoured ones with 
respect to the onset of the instability, typically having t^d ~ Id — 1 yr. On the other hand, however, too strongly magnetized 
and/or too rapidly rotating magnetars will spin down before the formation of their crusts, rendering this discussion moot. 

It is beyond the scope of this paper to provide a detailed analysis of the fate of the unstable, super-rotating core region. 
Nevertheless, we can speculate. Accepting that an initial magnetic equilibrium such as that in the left panel of Fig.[^is quickly 
rendered dynamically unstable by the super-rotation of the closed-field-lines region, the magnetic field will try to rearrange 
itself and find a new stable equilibrium. As long as the neutral line remains in the core, this does not seem viable; any new 
configuration will still contain closed field line regions which will again lead to a crust-core spin-lag. 

Thus, a stable equilibrium is likely to be reached only by the eviction of the neutral line and the attached toroidal field 
from the core and into the crust. We therefore suggest that the failure of the magnetic field to enforce complete crust-core 
corotation drives a magnetic field evolution in very young neutron stars akin to moving from the left to the right panel of 
Fig.[T] 

This hypothetical evolutionary scenario sees the crust as the natural depository of toroidal magnetic field and energy 
initially located in the core. As pointed out earlier, this kind of rotation-driven magnetic field evolution could be especially 
favoured in magnetars because these magnetized objects experience a rapid spindown soon after the formation of the crust. 
In other words, and according to the above suggestion, magnetars may have their strong interior toroidal fields pushed into 
their crusts in ~ 1 d — 1 yr after their birth. This prediction sits well with our theoretical understanding of magnetar thermal 
properties and seismic activity. A strong, evolving toroidal field in the crust can easily provide the energy required to heat 
magnetar surfaces to the observed levels ( Pons &: Geppert|2007 Kaminker et al.|2006 |Ho et al. |2012 (. The same process could 
be responsible for fracturing the crust and powering high energy events like magnetar flares ( Thompson fc DuncanpOOl (. 

At this stage it is more difficult to speculate on the repercussions of the proposed magnetic field evolution on canonical 
systems like radio pulsars, or the less magnetized central compact objects (CGOs). The smoother electromagnetic spindown of 
these objects suggests a much longer spin lag longevity, as t^d ~ 10® yr — 10"* yr. It should be recalled that the viscous coupling 
timescale can be of comparable length ( |Melatos]|2012[ ) and therefore corotation could be achieved regardless of magnetic 
coupling. It is more likely, however, that the Kelvin-Helmholtz instability will drive magnetic field evolution much earlier than 
that, especially if the internal toroidal field is comparable to the surface field, i.e. 
latter case, and accounting for superconductivity, the instability would set in after 


OV 

^0 


10®® G; see equation (67l. In this 


1 yr — 100 yr. 

Although the previous scenario is attractive for all the reasons discussed, it is worth reiterating that it is based mainly in 
conjecture. However, if for some reason it is not realised, it is possible that a persistent ‘quasi-steady’ turbulent state could be 
established, driven by shear flows at the interface between the two regions. Turbulent flows in the interiors of neutron stars 
have many potentially interesting observational consequences (for example see Mastrano fc Melatos||20()5 Melatos|2012 and 
references therein). Establishing whether the closed-field-line region is evicted from the core or not is a difficult task that we 
anticipate will require high-resolution numerical simulations to resolve. 

Our final remark concerns the role of neutron superfluidity in the evolution of the super-rotating core region. This could 
be relevant if the persistent spin lag survives until the bulk of the stellar core has become superfiuid (~ lyr — 100yr). Given 
the previous (very approximate!) instability timescales, this scenario is more likely to be realized in young pulsars, although 
the low magnetic field/spin frequency end of the magnetar distribution cannot be excluded. The neutron vortex lattice which 
enables the superfluid to rotate will generally thread both the core region coupled to the crust and the decoupled super- 
rotating region. The mutual friction force will efficiently couple the neutrons to the protons in both regions, giving rise to a 
differentially rotating superfluid. In order for this to be possible, the vortices that thread both regions must adapt their spatial 
distribution and, at the same time, bend in the vicinity of the boundary field line A (see Fig. [^. However, the self-tension 
of the vortices will prevent this from happening indefinitely. In fact, it might be possible that corotation is established in the 
entire core under the combined action of the magnetic field and the superfiuid vortices. This is an attractive scenario but is 
one that, as we have pointed out, requires an initially long-lived persistent crust-core spin lag. 


6 CONCLUDING DISCUSSION 

In this paper, we derive conditions for a magnetic field threading the interior of a neutron star to enforce long-term corotation 
between the slowly-spinning down crust of the star and the superfiuid core. We show that a magnetic field with open field 
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lines allows for complete corotation between the two components. However, somewhat counter-intuitively, the presence of any 
closed field-line region in the core of the star causes that region to be magnetically decoupled from the rest of the star. As the 
crust and the open-field-line region of the core spin down in unison, a velocity lag will build up between the closed-field-line 
region and the rest of the star. 

Example magnetic field topologies that will generate such persistent spin lags include the popular ‘twisted torus’ config¬ 
urations, in which toroidal field threads the closed-field-line region of the internal poloidal field (see the left hand panel of 
Fig. 0 . Our results indicate that, soon after the formation of the crust, a young neutron star with, for example, a twisted torus 
magnetic field will develop a persistent spin lag in the region enclosed by boundary A in the left panel of Fig.[^ The long-term 
fate of this super-rotating core region is unclear. The velocity difference between the two regions gives rise to a current, and 
a corresponding Lorentz force that needs to be considered when calculating new magnetohydrostatic equilibria. However, we 
expect any such equilibria to be unstable to e.g., the classical Kelvin-Helmholtz instability at the interface between the two 
regions (boundary A in the left-hand panel of Fig. [^. In Section]^ we speculate that any long-term evolution would involve 
the eviction of the toroidal component of the magnetic field into the crust, thereby allowing the entire core to corotate with 
the crust. 

Fully addressing the question of the long-term fate of the super-rotating core region is a challenging task that requires 
sophisticated numerical simulations that, we believe, are beyond the capabilities of current techniques. Some difficulties include: 
simultaneous resolution of the small-scale structure induced by the instabilities at the shear layer and the global magnetic 
field topology; and the range of timescales in the system, from the secular spin-down and Kelvin-Helmholtz timescales (~ 
days - years) down to the Alfven timescale governing the dynamical adjustments of the magnetic field. We therefore advocate 
a multipronged approach, include calculations of quasi-static, magnetohydrostatic evolutions to understand the influence of 
the poloidal Lorentz force on the toroidal-field-line region, as well as full magnetohydrodynamic simulations to understand 
the relevant instabilities at the shear layer. 

Although our analysis makes no particular assumption about the properties of neutron star matter (equation of state, 
degree of stratification etc.), it does involve two main simplifications with respect to the magnetic field. The first one is the 
assumption of a dipolar structure for the poloidal field which facilitated the full use of the magnetic coordinates formalism. 
Of course, these coordinates remain a well-defined notion and can be employed for the description of a general axisymmetric 
field, but this may come at the price of losing a large part of the analytic elegance of the dipole case. Although we have 
not pursued this issue further, it is unlikely that any of our conclusions will change for a magnetic field of more general, 
axisymmetric structure. This may not be true for non-axisymmetric fields where, for example, the poloidal field is misaligned 
with the stellar rotation axis | Lasky fc Melatos|2013 1, however we suspect any field topology that admits closed field lines in 
the core will exhibit the same phenomenology. The second simplification is the use of the ‘classical’ Lorentz magnetic force 
in the Euler equation. Strictly speaking this is ‘wrong’ given that the bulk of the proton fluid is expected to condense to a 
superconducting state around the same time as the formation of the crust. It is well known that the magnetohydrodynamics of 


superconducting neutron stars can be quite different to that of their ordinary matter counterparts (Mendell 1998 Glampedakis 
et al.|20Tl I. One of the key difference^ is the form of the magnetic force, which in a superconducting system originates from 


the smooth-averaged self-tension of the quantized protonic fluxtubes. Future work should address this approximation of our 
model by using the full machinery of superconducting magnetohydrodynamics. 
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APPENDIX A: CRUST-CORE COUPLING FOR A NON-UNIFORM DIPOLE FIELD. 


Here we consider the non-uniform dipole field /(r) = C combined with the assumption of uniform density. Then, eqn. (281 
can be solved with separation of variables, leading to the solution 

f1T/2 


y^{e) = - —= ( -Ppfl -t- 2Z3pnDnDnp ) / M' Vsin . 

AVsind ' ' Je 


Using this result and Bp = 2AC cos 6 in the formula for the torque, ALc, we have after some rearrangement: 

Ncc = —7pD -(- 2ZIDnDnpirL, 


(Al) 

(A2) 


^ Another difference is the “boosted” superconducting Alfven speed = HcB/A.ivpp-, this expression has been used in some of our 
estimates in this paper. 
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with the effective moments of inertia 


rTr/2 


.-.2 ^ COS 


7x = 27ri?cpx / d9 sin 6—^ _ 

VsinS 


a r '^/2 

— I dOWsine'. 


Jo 

On the other hand, the total moment of inertia for each fluid is 


le 


/■ 


7x = / dV sin^ 6 = An 


n -K j 2 


drdOr'^ sin^ 9p^. 


For uniform density this simplifies to 


7x = ^7?®Px. 
15 


(A3) 


(A4) 

(A5) 


Our final task is to check the value of the ratio 7x/7x. We have, 

d9 sin®^^ 9 cos 9 

Direct numerical integration reveals that to numerical precision I^/ 
electrons) with the crust. 



r/2 


d9Wsin9'. 


(A6) 


0 

= 1 which implies a fully corotating core (protons and 


APPENDIX B: MAGNETIC COORDINATES 


For a given axisymmetric poloidal field, Bp = x Vp, we can set up a system of magnetic coordinates, {■)/), x, p}, where x 
measures distance along a field line (without loss of generality we can set x = 0 at the magnetic equator). 

The resulting coordinate system is orthogonal but curvilinear with line element, 

dr ■ dv — d'lp'^ + dx^ + (Bl) 

where are all functions of and x, and it is easy to see that h^p = w. By construction. Bp is tangent to ip surfaces, 

i.e.. 


Bp = Boiip,x)x- 

The grad operator takes the form (where as always a ‘hat’ denotes a unit vector) 


II'iIj '^X fl'ip 


We then find, 


Bp = Vi/> X V<p = 


-X- 


(B2) 


(B3) 


(B4) 


hiphp 

For the remainder of the analysis we assume a dipole field, implying the stream function can be expressed in the general 


form 


i,{r,9) = Af{f)sin^9, (B5) 

where f = r/R. As we have seen, the components of this field in spherical coordinates are 

BI = 2aJ- cos 9, B^ = -aJ- sin 9, (B6) 

Ttr 

where /' = df{x)/dx. 

From this we can obtain the metric component h^: 

hp (r, ^ g [ 4/^ cos^ 9 + (rf'f sin^ 9 ] . (B7) 

Expressing the line element for the two-sphere in terms of h,p and h^, we can obtain three equations 

1 “ (^/O^sin^^i [4/^cos^6»-|-(r/')^sin^6l] (B8) 


(^x^) = r'^ - 4: f cos^ 9 [4 f cos^ 9 + {rf'f sin^ 9] ^ j , (B9) 

=—2sm9 cos9rfff'[4f^ cos^ 9 + {rf'f sin^ 9] (BIO) 

It is straightforward to see that the bottom equation is trivially satisfied by the first two. Simplifying the first two 
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equations and taking square roots, 

dX _ 1 r,j-COs6» ^ j,2 _ 2 „;„2 fli-i/a 




X 

dX _ , -./Sing 
89 ^ K 


[4/ cos^6l + (f/') sin e]’ 


[4/^cos^6l + (f/')%in^6l] 


- 1/2 


(Bll) 


(B12) 


Somewhat trivially, drX > 0 for all for n dipole poloidal field, and we therefore consider only the ‘-I-’ sign in equation 

(Bill. On the other hand, dex oan be positive and negative along any given field line. For example, consider field line A in 
the left hand panel of Fig. As one moves from the innermost equatorial point on A in the clockwise direction, dex > 0 
until the apex of the field line is reached, at which point dex changes sign, and becomes negative as one moves further in the 
clockwise direction. We accordingly introduce e = ±1, 0 into equation ( |B12[ ). 

Motivated by solutions in the case when /(f) = f“ (see below), we search for a coordinate transformation of the form 

x{r,0) = G{r)cos6. (B13) 

Inserting this ansatz into equations (|B11[| and (|B12[| gives respectively 


dG 2/ 


dr h 


= -A [4f cos^ 9 -f (rf'f sin^ 9 


- 1/2 


(B14) 


G = —erfj— [4/^ cos^ 9 + sin^ 9] 


(B15) 


Dividing equation ( |B14[ ) by (B15I eliminates the dependence on 9, and gives us a first-order ordinary differential equation for 
G. The formal solution of this equation is 


G = exp(-?/lfi) 


which can be put back into equation (B13l, leaving the final coordinate transformation: 

cos 0. 


X(r,0) =exp(^-^ f 

Inserting this back into equation (|B8[) for the metric coefficient, h^, yields 


hxir,S) = -^rrf' exp J df^-p^ [4/^ cos^ 0-|-(f/')^ sin^ e] 


- 1/2 


(B16) 


(BIT) 


(B18) 


For a given dipole field, one prescribes /(r), and can derive the full coordinate transformation between spherical polar 
coordinates and magnetic coordinates using equations (B5l and (BlTl. 

As an example of this we consider a simple power law, /(f) = f“, which includes a uniform field (a = 2) and the ‘standard’ 
dipole field (a = 4) as special cases. In this case, the coordinate transformation becomes 

'tl){r,9) = Af^-siC 9, x{r,9) = Gf-^^°'"cos9, (B19) 

where C is a constant of integration. 


APPENDIX C: SOME GEOMETRIC APPLICATIONS OF MAGNETIC COORDINATES 


As a first application of our magnetic coordinates, we calculate the infinitesimal surface element dS on a sphere of radius rg. 
This is given by 

dS = h^dtp [h^dtl)^ + h^dx^Y^^ ■ (Cl) 

On the spherical surface we have dx = —{d^ro/d^ro)d'ip and then 


dS — hiph'ijj 


1 + 


h^d^r 


1/2 


dipdij). 


(C2) 


As a second application of the magnetic coordinates, we derive the equation for a circle of radius ro = ro(x, i/). Beginning 


with the stream function, 'ip = Af{f) sin"^ 9, we use equation (BITI to eliminate the sm 9 term: 


A 


G2 


(C3) 


which is an implicit equation for the circle since G = G(ro) and / = /{tq/R). From this we obtain the partial derivatives 


drp 

dip 


1 

A 


r 


R 


X 


G2 


:G' 


^ +2/x^7^ 


G3 


drp ^ yjc 
dx G2 


/' 


R 


X 


G2 


G' 


^ +2/x^ — 


G3 


(C4) 
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Their ratio is 


and we also have 


From these results we find 


d^ro _ G 
dx^o 2Afcos6’ 


hij] 


-A sin 61 


erp/' 

G 


hxd^ro _ e //' 
dxTo 2 / 


In order to make contact with the main body of the paper we also need to calculate the inner product 


f •x = ^ = 2/cos6»[4/cos^6i +(f/')^sin^6»]"^/^ 
B-p 


We also have 


f ■ %j) = ±[1 — (f • = iff' sin0[4/^ cos^ 6 + {ff'Y sin^ 0] 




It is then straightforward to verify the identities: 


(r-x) 


1 + 


hxd^rp V 
h^dxTo ) 


11/2 


= 1 , 


[r-Y) 


1 + 


hxd^rp Y 
h^dxTo ) 


11/2 


1 f f' 
-±-tan6>^. 
^ / 


For the example given in Appendix]^ f = , the above formulae reduce to 

d^ro Ra /ro\-“-2/“'= (hxd^rpY^ - o q\a 2 ^, 2 . 2 ai-i /2 

— -n , , „ — = — tan 6, r ■ x = 2 cos 9 4 cos 9 + a sm 6* 

V if / \ h^OxTo / 4 '■ ■' 


dxTp 


2 AG cos 9 


and from these the identity (CIOI follows directly. 


(C5) 

(C6) 

(C7) 

(C8) 

(C9) 

(CIO) 

(Cll) 


APPENDIX D: THE POLOIDAL EULER EQUATIONS 


The poloidal component of the perturbed Euler equation ([^ is (after omitting the Coriolis term and the unimportant ~ B' 
mutual friction term): 


p'i!p = —Vp (Bo • b) + [ (b • V)Bo + (Bo • V)b]p . 

Using magnetic coordinates, we have for the individual terms of (IDlIl 


[(b.V)Bo]p=i)> 
[(Bo-V)b]p =!/> 


h h , , Otjjhip 

rixjjrix nipitip 


1 

r B^b^ ^ 1 

+ x 



Bpby 


and 


Bo • Vbp - Y^dphx - 

.'iprix riipiL^p 


Ip 


+ x 


„ , Bpbp ^ , B‘^b‘^ ^ , 

Bo • Vbx + -l—r-ophx - -j—r-dxK 

iL'ijjIl'x I Lx* 


dpP \ , X 


Vp (Bo • b) = ( B'^b'^ + ^ J + ^dx ( B^b'^ + 


dpP 

h2 h2 I ■ 


Recalling the expressions, 

T = mb‘^, B°p=IV^, hp^VpxV^, 

then after some straightforward manipulations we arrive at the following components of ( |D1[ ) 

^ ^ dlP-KxdxP + 2Tx’^d^hp + 2^dphpdpl3] 


-4'Kppdp'i'p = dp ['^ + 


hi ^ hlhl j ^ {hxh^Y 




A nr I T"! dpP 

-47rppdx4'p = dx[ -rpr + 


hl hlhl ) (hph^Y 




dxdi^P + K^dxP - 2TX-^dxK 1 , 


where 


^X^X ^X^i/’ 


Kx = ^ + 


hip hx hip 


Kit = 


diph(p Oiphip diphx 


(Dl) 

(D2) 

(D3) 

(D4) 

(D5) 

(D6) 

(D7) 

(D8) 
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